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ABSTRACT 

Aims. The On-The-Fly (OTP) imaging technique enables single-dish radio telescopes to construct images of small areas of the sky 
with greater efficiency and accuracy. 

Methods. This paper describes the practical application of the OTP imaging technique. By way of example the implementation of the 
OTP imaging technique at the NRAO 12 Meter Telescope is described. 

Results. Specific requirements for data sampling, image formation, and Doppler correction are discussed. 
Key words, telescopes ~ methods: observational - techniques: image processing 
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1. Introduction 

In an elfort to become scientifically more efficient, astronomi- 
cal observatories have incorporated innovative observing tech- 

Q niques to increase the throughput of their telescopes. The On- 
The-Fly (OTP) imaging technique is an innovation which has 
been adopted at a number of radio observatories over the past 
^ ^ foity years in a variety of forms. In the 1960's and 1970's, ra- 
dio observatories which operate at centimeter wavelengths im- 
plemented the "drift scanning" technique (sometimes referred to 

J> as "nodding") to acquire continuum imaging measurements (cf. 
|Haslam et a/.| ( |197"0 )). Drift scanned measurements are obtained 
by slowly slewing the telescope in elevation with the Earth's ro- 
tation providing the second dimension. In the 1980's, a number 
of radio telescopes added beam-switched continuum imaging ca- 
pabilities. Beam-switched continuum observations are acquired 
by slewing the telescope in a two-dimensional raster pattern over 
r — a given region of sky while rapidly wobbling the telescope's sub- 
reflector to obtain differential measurements of the total power. 

^ The [Emerson et ar| ( |1979| l [EKH] algorithm allowed the decon- 
• ^ volution of these beam-switched measurements into total power 
images. This powerful imaging capability added considerably to 

J-H our ability to understand the centimeter- and millimeter-wave 

^ continuum structure of the interstellar medium. In the 1990s a 
number of radio observatories implemented both spectral line 
and continuum OTP imaging capabilities {cf. Mangum (jT999]l 
and |Mangum et a/.| ( |2000|| l). 

As with many innovations, this "new" OTP imaging tech- 
nique was just an adaptation of an old observing technique. The 
observations which signaled the birth of radio astronomy, made 
by Karl Jansky in 1932, were actually OTP observations (see 
|Janksy| ( |1932| l and Figure [TJ. Jansky's 20.5 MHz synchrotron 
measurements of the galactic center represent the first astronom- 
ical OTP observations. In the following, we describe the adap- 
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Fig. 1. Jansky and his telescope ( |Sullivan| ( |1978| )). 

tation of Jansky's imaging technique to total power images ob- 
tained using a single dish telescope. 

2. Advantages of OTF Observing 

OTP mapping is an observing technique in which the telescope 
is driven smoothly and rapidly across a region of sky, or "field", 
while data and antenna position information are recorded con- 
tinuously. This technique is in contrast to traditional mapping of 
discrete positions on the sky, which is sometimes called "step- 
and-integrate" or "point-and-shoot" mapping. The advantages to 
OTP mapping are: 

1 . Telescope overhead is reduced significantly, since a specific 
position on the sky does not have to be acquired within a 
given tolerance (i.e. minimizing "dead-time"). 

2. The entire field is covered rapidly, minimizing changes in the 
properties of the atmosphere and the system, including an- 
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tenna pointing and calibration. Systematic changes may oc- 
cur from map to map, but such effects average down rapidly 
and may be correctable by cross correlation techniques. In 
general, global changes from map to map are more benign 
and easier to correct than drifts across a single map field. 
Higher observing efficiency: The radiometer equation for 
switched measurements is given by the following 
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where cr is the RMS noise in a total power measurement, 
Tsys is the system noise temperature, Av is the spectral reso- 
lution of the measurement, and rjspec is the spectrometer ef- 
ficiency. The optimum duration of an OFF measurement for 
any switched measurement is given by { fial\\l916) ) 
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where N is the number of ON measurements made per OFF 
measurement, such that tscan - toff + Nton. Therefore, for 
optimal observing the RMS noise in a measurement is given 
by 



cr - 



T 
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where tscan, as defined above, is the total ON plus OFF mea- 
surement time. For N>100, cr is approximately a factor of 
two smaller than it is for an equivalent position switched 
measurement, which leads to a factor of four improvement 
in observing time efficiency. 

3. Sampling 

3. 1 . Detection and Sampling: The Spectrum of tiie Detected 
Noise and Signal 

When conducting single dish imaging observations, it is impor- 
tant to keep in mind the following facts about sampling and alias- 
ing in radio astronomical mapping data. If you want to represent 
the full resolution of the telescope, you have to sample the data 
often enough to represent all the spatial frequencies detected by 
the antenna. Figure |2] illustrates a typical continuum detection 
scheme of a single dish telescope; this is the simplest case, but 
the same principles apply to detection in a spectrometer, or even 
direct digitization at the IF and subsequent processing by soft- 
ware or firmware. In the example shown, the radio astronomi- 
cal signal is amplified and bandpass filtered, before passing into 
a square law detector, possibly some simple filtering, and then 
into an Analog-to-Digital (A/D) converter The A/D converter, 
or later processing, usually incorporates some signal integration 
before delivering the sampled data. 

Figure |3] illustrates the low frequency spectrum of the sig- 
nal and noise immediately after the square law detector and in 
the subsequent processing steps. With the antenna scanning at a 
constant speed over a point source, its response produces, after 
the detector, a time-varying signal whose spectrum corresponds 
to the spatial frequency response of the antenna. The horizon- 
tal axis of Figure |3]could be spatial frequency, measured in units 
proportional to ^ (d is a baseline distance), or the frequency after 
the detector measured in Hz. Any antenna has a spatial frequency 
cut-off; with an antenna of diameter D and scanning speed of O 
radians/second, after the detector the signal voltage will have a 



Fig. 2. Typical continuum detection scheme for a single dish 
telescope. 

Single Dish Spfitifil Frequency Response 




Relative spatial frequency 

Spatial frequency response of antenna 

Spatial frequency cutoff of antenna 

Detector output spectium before integration 

Possible response of simple RC filter after detector 
Square bos integration, Nyquist sampled 

Fig. 3. Illustration of the low-frequency spectral response of the 
signal and noise after square law detection. 

cut-off frequency of (j)^^, now measured in Hz. This is illus- 
trated by the bold curve in Figure[3] The signal cut-off frequency 
is shown by the vertical dashed line. 

Figure |3] also shows the spectrum of noise (the dashed line) 
immediately after the detector; this is strictly the autocorrelation 
function of the IF passband, which means in practice that within 
the narrow range of post-detection frequencies of interest, it is 
white noise, with a flat spectrum. It is normal to include a sim- 
ple low pass filter immediately after the detector, shown sym- 
bolically by the RC time constant in Figure |2] This might have a 
gentle frequency cut-off such as illustrated by the dash-dot line 
of Figure |3] 

The digitizer usually includes some time integration as part 
of its operation; the output samples might be sampled S times per 
second, but each sample is an integration of ^ seconds of data. 
This integration is equivalent to convolving the time sequence of 
data with a square-box function - equivalent also to convolving 
the sky image along the scanning direction by the corresponding 
square-box function. The integration of S seconds by this square- 
box translates into an average of 5 * Q radians of angle along 
the sky scanning track. The time-domain spectral response, or 
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equivalently the spatial frequency response, of this square-box 
integration is illustrated in Figure |3] for the case where the sig- 
nal is exactly Nyquist sampled (two independent samples per ^ 
interval on the sky). 

The ideal signal processing and detection scheme might use 
a perfect square-edged anti-aliasing filter in place of the sim- 
ple RC-filter illustrated in Figure |2] followed by a sampler op- 
erating at least at twice the spatial frequency cut-off, with an 
infinitely narrow window in the sampler, thereby avoiding the 
loss of high frequencies caused by the normal square box inte- 
gration. However, to avoid the complexities (and impossibility) 
of building a perfect anti-aliasing filter, very often in practice the 
square-box integration is retained, but with a substantial increase 
in the sample rate above the Nyquist value. 

3.2. The Consequences of Undersampling 

The spectrum of noise and of signal at the output of the detector 
of a total power system depends on the details of the electronic 
design, but a representative system is illustrated in Figures|2]and 

s 

With OTF observing, the astronomical sky is convolved with 
the telescope beamshape. The telescope acts also as a low pass 
filter, removing all astronomical components beyond the intrin- 
sic angular resolution of the antenna. The telescope beam and 
telescope scanning speed together give a low pass filter, with the 
conversion from spatial frequency to temporal frequencies being 
given by 

R(f) = r(sf) * Q, (4) 

with R(f) being the temporal frequency response in Hz and r(sf) 
the spatial frequency response measured in units of j, with O 
the telescope scanning speed in radians/second. The respective 
responses are illustrated in Figure [3] 

3.2.1. Aliasing 

The minimum sampling interval on the sky, for an antenna with 
maximum dimension D, is at intervals of Assume that one 
undersamples on the sky, rather than later in the data processing. 
Suppose you have a 10m dish, but you only sample at /l/(2x8 m) 
rather than the A/ (2 x 10 m) that you should. This means that 
the spatial frequencies present from the dish baselines of 8 m to 
10 m get reflected back into the spatial frequencies of 8 m down 
to 6 m. Not only have spatial frequencies from the 8 m to 10m 
baselines been lost, but valid spatial frequencies from baselines 
of 6 m to 8 m have been coiTupted. You can't tell if structure 
in your map with a spatial wavelength of A/1 m is genuine, or 
was really structure at A/9 m which has been aliased on top of 
any genuine A/7 m spatial wavelength signal. In this sense, un- 
dersampling the sky is really twice as bad as you might have 
thought. 

Figures |4] and |5] show the spatial frequency response and 
aliased noise power for Nyquist sampling while using a square 
box sampling function, of width equal to the sampling interval. 
The power Pj in the random noise fluctuations in frequencies up 
to the properly sampled spatial frequency, indicated by the verti- 
cal dashed line at 0.5 on the abscissa of Figure[4j is proportional 
to: 




Spatial Frequency Response 




Spatial Frequency (1/t) 

Fig. 4. Spatial frequency response for square box Nyquist sam- 
pling. 



Aliased Frequency Response for Nyquist Sampling 
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Fig. 5. Aliased noise power for square box Nyquist sampling in 
the absence of a separate anti-aliasing filter. 

The power A„ in the random fluctuations in frequencies be- 
yond that frequency, in absence of any anti-aliasing or other fil- 
ters, is proportional to: 

A„=r(^)^. (6) 

Jo.5 \ j 

This noise power A„, which cannot contain any astronomical 
information, is then aliased back into the signal passband below 
0.5. The random (RMS) noise voltage is then increased by the 
ratio: 

^^^ = 1.137 (7) 

equivalent to a loss in signal-to-noise of 12%. Therefore, the 
noise at frequencies to the right of 0.5 on the abscissa of 
Figure |5] aliased back into the region of interest to the left of 
0.5, leads to an overall degradation in signal-to-noise of the im- 
age of 12%. 

The noise power degradation due to sampling is a function 
of the level of oversampling in the image plane. Figure |6] shows 
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Signal— to— Noise Ratio and Sampling 
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Fig. 6. Signal-to-noise degradation due to noise aliasing of an 
inadequately sampled image using square box sampling. 



the amount of signal-to-noise degradation one can expect due to 
noise aliasing of an inadequately sampled signal using square 
box sampling. 

How important undersampling is depends on exactly what 
the illumination taper is, how important it is that you retain the 
maximum possible resolution of the telescope, how good a dy- 
namic range you want in the observations, and at some level how 
much fine scale structure there is in the source itself. If you only 
sample at 0.833 Nyquist (e.g. FWHM/2 rather than FWHM/2A 
for a -13 dB illumination taper), what matters is the energy in 
the data at spatial wavelengths shorter than /l/(2 x 0.8 x D). So in 
a sense you need to ask what the illumination taper is at a radius 
of 0.4 X D on the dish surface. The spatial frequency response 
of a single dish is the autocorrelation function of the voltage 
illumination pattern. So, you need to calculate how much area 
there is under the 2-D autocon^elation function beyond spatial 
frequencies of 0.8 x D, compared to the area within 0.8 x D. This 
ratio is some measure of the dynamic range. A better definition 
of dynamic range might take into account the spatial frequency 
structure of the source. If the source has no structure on scales 
smaller than A/ (2 x 0.8 x D), then you don't need to sample at 
the full A/ (2 X D) anyway. 

There are circumstances where it is perfectly rigorous to un- 
dersample the data in the absence of a separate anti-aliasing fil- 
ter For example, if you have a 10 m dish, and you are taking 
data to compare with other observations that were made using 
aim dish at the same wavelength (or the equivalent number 
of wavelengths at some other frequency) then you only need to 
sample the data at A/(2 x 5.5 m) or /l/ll m. This is so because, 
when sampling a 10 m dish as if it were a 5.5 m dish, the spa- 
tial frequency components from baselines of 5.5 m out to 10 m 
will be reflected back into the data as if con^esponding to base- 
lines of 5.5 m down to 1 m. So, the spatial frequency terms of 
the 1 m baseline and below will not have been corrupted. The 
data analysis of this undersampled data would apply a spatial 
frequency cutoff at 1 m, and there will have been no corruption 
in this smoothed data caused by the undersampling. Putting it in 
more general terms, if you are going to be smoothing observa- 
tions from a dish of diameter D to simulate observations made 
with a smaller dish of diameter d, then the sampling interval only 
needs to be A/(d + D) rather than Al{2 x D). 



3.2.2. Beam Broadening 

Since the telescope beam is effectively convolved with a square- 
box function equal to the distance the telescope moves in one 
sample period, there is a degradation of spatial resolution due to 
the data acquisition process. If you set an arbitrary criterion that 
the telescope beam shouldn't be broadened by more than 1%, 
then you need to sample at least 4.5 points per FWHM beam 
width. Note that since Nyquist sampling is about 2.4 points per 
FWHM, assuming an illumination taper of about -13 dB; for 
a uniformly illuminated circular aperture the factor would be 
about 2.06, while for a -20 dB taper it becomes 2.60. This means 
that one needs to sample at about twice Nyquist to avoid degra- 
dation of resolution worse than 1%. 

3.2.3. Sampling Summary 

1 . To avoid both beam degradation and loss of signal/noise ra- 
tio, one must sample at least twice Nyquist, or about 5 points 
per FWHM. 

2. With few exceptions, undersampling at data acquisition is 
not recommended. 

3. There are other aspects that make it desirable to sample more 
often than the Nyquist rate, as is often recommended for 
OTF observing at many millimeter-wave telescopes. These 
are practical points like how well gridding or interpolation 
works with a finite sized gridding or interpolation function. 
For example, a little oversampling may enable you to reduce 
the convolution (interpolation) function by a factor of a few, 
saving a huge amount of computational overhead at the ex- 
pense of a few per cent more data. 

4. Gridding (Convolution) 

The final image is constructed by smoothing the data at whatever 
coordinates they were observed and then re-sampling them on a 
regular image grid. The smoothing is, in fact, an interpolation, 
which can include weighting and averaging, rather than a strict 
convolution. However, when the data are sampled in a dense and 
nearly uniform fashion, the interpolation approximates a convo- 
lution. Figure|7]shows some representative convolving functions 
which may be used in the data gridding. The eff'ect of the convo- 
lution is best seen by examining its impact, a multiplication, in 
the Fourier, spatial frequency, space. The Fourier transforms of 
the convolving functions are illustrated in Figure |8] 

The spatial frequency response pattern of the single dish 
is eff'ectively multiplied by the Fourier transform of the con- 
volving function. If the function were to have no effect on the 
data, then its Fourier transform should have value 1.0 out to 
the maximum spatial frequency of the telescope and 0.0 outside 
that radius. This function is, in one dimension, a sine function, 
sm(nx/NS)/{7Tx/NS) and, in two dimensions, a Bessel function 
Ji{7Tx/NS)/{7Tx/NS). Unfortunately, both of these functions die 
off very slowly, making the convolution operation excessively 
expensive. A compromise function is obtained by tapering the 
convolution function by a fairly narrow Gaussian, allowing the 
function to be truncated outside a restricted region. In the Fourier 
space, this has the effect of reducing the response to the largest 
spatial frequencies available to the telescope. These spatial fre- 
quencies are probably already tapered by the response pattern of 
the feed horn, which is usually designed to avoid sensitivity to 
ground spillover 

Since we have a circular telescope and are gridding 
the image in two dimensions simultaneously, we use the 
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Fig. 7. Representative convolution functions. 



Fig. 9. An example of a raster scanning pattern. In this example, 
three scanning rows are acquired per reference ("off") measure- 
ment. 
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-0.2 0.2 

Image Diameters 

Fig. 8. Representative Fourier transforms of convolution func- 
tions. 



circularly symmetric version of the convolving function, 
Ji{r/a)/(r/a)exp(-(r/by). Optimizations of this function give 
the default values for a, b, and c of a = 1.55 (single-dish-beam- 
size / 3), b = 2.52 (single-dish-beam-size / 3), c = 2 (see Table[T i, 
and a support radius equal to the single-dish-beam-size (FWHM) 
( |Schwab| ( fT98 0^). 

It is important to choose the sampling of the output image, 
and the width and support of the convolving function to avoid 
aliasing. This topic was discussed at length in ^3.2.1 A suitable 
cell size would be less than or equal to one-third of the single- 
dish-beam-size. 



Other convolving functions have various undesirable ef- 
fects. For example, a sinc^ has a triangular Fourier trans- 
form. All non-zero spatial frequencies are weighted down 
linearly with frequency, degrading the telescope resolu- 
tion. All other positive-only functions have a similar ef- 
fect. One can construct convolving functions such as 
sin(.x/a) sin(y/a)/(x/a)/(y/a)exp(-((x/b)^ + (y/b)^)) which will 
actually improve the spatial resolution of the image over that of 
the single-dish beam. They do this at the cost of down-weighting 
low spatial frequencies and thereby reducing the signal-to-noise 
of the observation ( |Greisen| ( |1998| l). 




Fig. 10. An example of the spiral scanning pattern. 



5. Scanning Geometry 

The scanning geometry used to acquire an OTF map is usu- 
ally dependent upon the type of source you are imaging and 
the ability to command the antenna to execute complex patterns. 
Scanning geometries that have been implemented at existing ra- 
dio telescopes include: 

Raster: Usually the first scanning geometry implemented at a 
telescope as it is the simplest to implement within a telescope 
monitor and control system. Figure|9]is a pictorial display of 
this scanning geometry. 

Spiral: An efficient scanning geometry when imaging a circular 
or elliptical source (see Figure 10 1. 

Hypocycloid: For imaging a non-circular source which is large 
in relation to the resolution element of your observations, 
the hypocycloid is a very motion-efficient pattern (i.e. no 
antenna motion efficiency losses due to turn-around; see 
Figure 111. 



6. OTF Imaging in Practice 

It is often important to estimate the amount of integration time 
and the RMS noise level an OTF image will have based on the 
observing parameters. In the following, we derive these quanti- 
ties for a raster scanned OTF map. We also discuss an often over- 



J. G. Mangum, D. T. Emerson, and E. W. Greisen: The On The Fly Imaging Technique 




OTF Row Sampling Rate 



Fig. 11. An example of the hypocycloid scanning pattern. 

looked aspect of wide-field spectroscopy involving the doppler 
correction. 



6.1. OTF Map Parameters 

In the following we describe some of the OTF map parameters 
that one must generally calculate in order to setup an OTF ob- 
seiTation. 



6.1.1. Row Sampling Rate 

It is very important that you set up your map to be properly sam- 
pled in all coordinates. If you under sample, you will miss infor- 
mation in the image field, you will be unable to combine your 
image accurately with that from interferometers or other single- 
dish telescopes, and you may introduce artifacts from the anal- 
Keep in mind that you can always 



ysis algorithms (see 5 3.2.1 



smooth the map after it is taken to degrade the resolution and 
improve signal-to-noise. 

The scanning rows must be spaced at no more than the 
Nyquist spacing, which is also known as "critical sampling" and 
is given by 



'lobs 

~2D 



radians 
30918 



D{m)votAGHz) 



(8) 



In practice, though, a small amount of oversampling is recom- 
mended. If the rows are critically spaced, small scanning errors 
can result in the map being under sampled. In addition, the tail 
of the gridding function used by the analysis procedure extends 
slightly beyond the information cutoff" of the telescope, which 
will result in some noise being aliased into each grid cell (see 
5 3.2. 1 1. A reasonable approach is to sample both sky coordinates 



by: 



\2DI 




frequency (GHz) 

Fig. 12. Orow as a function of observing frequency for a 12m aper- 
ture. 
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D{m)Vobs{GHz)] 



- 6 arcseconds 



(9) 



where the factor of 0.9 is an oversampling factor and the factor 
of 5 is a guard band to accommodate any scanning errors. For 
the NRAO 12M Telescope 5 — 2 was used. Figure 12 is a plot of 
this relationship assuming a 12m aperture. 

6.1.2. Scanning Rate 

Map sampling along the scanning direction is dependent upon 
the integration (or "dump") time of the data acquisition sys- 
tem {i.e. spectrometer) being used. Since the spectrometers must 
integrate for a finite interval before being dumped, the data 
are "square-box" binned. To avoid noise aliasing and beam- 
smearing problems due to this binning (see §3.2.1 1, one needs 



to oversample in the scanning direction. The scan rate may be 
written as 



R(arcsec/ sec) 



(10) 



ump 



X 0.9 - (5 



where 6m is the Nyquist spacing (Equation[8]l, is an oversam- 
pling factor, and t^ump is the spectrometer data dumping inter- 
val. At the NRAO 12 Meter Telescope, fj,„„p =0.1 seconds for 
spectral line OTF, 0.25 seconds for continuum OTF. The general 
recommendation to observers when using the 12m system was 
to use a minimum value for Hos of 2, which will result in < 3% 
increased noise as a result of aliasing and minimal beam broad- 
ening (see 5 3.2.2 1. Note that it is quite acceptable to oversample 
by much larger factors than the minimum value suggested above, 
particularly in the scanning direction. There are some tradeoff's 
to be considered, however: 

1 . Large factors of oversampling will produce better signal-to- 
noise on a single coverage of the field and may result in sim- 
pler data processing, i.e. fewer maps to average together 

2. On the other hand, the longer a single coverage takes, the 
more susceptible your image will be to drifts in the system, 
such as pointing and atmosphere. This eliminates some of 
the advantages of fast mapping mentioned above. 

A good compromise is to use a row spacing computed from 
Equation |9] and a scan rate calculated from Equation [TO] with 
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rios > 2. Other circumstances may come into play, of course. For 
example, if the field size is very small, you may choose to slow 
the scan rate (increase rios) substantially. 

6.1.3. Map RMS 

When planning an OTF observing run you must calculate how 
much integration time you need for each sampling cell to reach 
the desired signal-to-noise level. First determine the integration 
time required using the standard Radiometer Equation (Equation 





O'cell - 



sys 



'{spec 



1 -I- 



tcell 
'off 



(11) 



where Tsys is the system temperature, risp^c is the spectrometer 
efficiency (typically ^ 1), Av is the spectral resolution, tcdi is the 
time spent integrating on a cell, and is the off-source inte- 
gration time. For spectral line OTF, the time spent integrating on 
a given cell is typically one second or less, while the off-source 
integration time is often several seconds (10 seconds was the de- 
fault value for the NRAO 12m Telescope). Therefore, tceii <s: toff 
and we can rewrite Equation 1 1 as 



1 /r,,, 

tcell - -7- 

Av \cr„// 



(12) 



The integration time in a map cell is a function of the scan- 
ning rate and the number of coverages of the map field, and 
should be calculated with respect to the time accumulated in a 
Nyquist sampling cell. For illustration we assume a raster scan- 
ning geometry. If we define 0^ to be the angular distance along 
the scanning row, 0iamp to be a "ramp-up" distance at the begin- 
ning and end of each row, and R to be the scanning rate (see 
^6.1.2[ ), then the time to scan one row is 



0JC + 20rc 

R~ 



If we also define 0, to be the angular height of the map and 0,-ow 
to be the angular sampling interval between scanning rows, then 
the number of rows in the map is given by 



The number of independent sampling cells in the map is 

(Ox + 201-amp) X 0y 



N, 



cell 



01 



(14) 



(15) 



where we are approximating the cells as rectangular. The inte- 
gration time per cell, is related to the above three quantities 
as follows 



trow^ row 

tcell = — X T] 

^cell 



(16) 



where 77 is a convolution function correction factor which takes 
account of the fact that we are using a tapered function (a 
Gaussian-tapered Bessel function) instead of a pure sine. We de- 
fine the convolution function C(u)„ as follows: 



C{u)m - I C{u)du 
Jo 



(17) 



The I] correction factor is then given by 



C{u\ 
C(0)„ 



(18) 



for the truncated functions used in the AIPS task SDGRD (see 
^and Table [1). Combining Equations [13] [14] [B] [16] and [18] 
we find that 



tcell - 



3.030^ 

R0rnw 



(19) 



Combining this relation for tceii with our Radiometer Equation 
we find that the RMS noise per cell in an OTF 
by 



(Equation 11 



map IS given 



CTcell 



R0n 



3.03Av/ 



(20) 



Note that a given RMS noise per cell obtainable with multiple 
coverages and/or polarizations is given by 



O'cell 



(21) 



where A^^. is the number of coverages and/or polarizations you 
combine to produce your final image. Note that you can also 
improve your map RMS by spatially smoothing it after gridding 
(see 



6.1.4. Total Map Time 

Again assuming a raster scanning geometry, the total time re- 
quired to acquire a map must include not only the integration 
time scanning across the field but also the time for cahbration 
and OFF integrations. This can be written as 



(13) t,o,^Nn 



toff 



teal 



(22) 



where Nrow is the number of rows in the map, trow is given by 



Equation 13 above, toff is the OFF integration time, t^ai is the 
calibration mtegration time (i.e. amplitude calibration and sky 
sample times), Nrpo is the number of rows per OFF measurement 
you do. Nope is the number of OFF measurements per calibration 
measurement (amplitude calibration and sky) you make, and e is 
the "overhead" time, which is the time that the telescope spends 
doing things other than integrating (like moving from OFF posi- 
tions to the map field). The value for e often depends mainly on 
how far away the OFF position is from the map, but 10 seconds 
is probably a good average estimate for vintage- 1980 millime- 
ter radio telescopes. Most of the time the amplitude calibration 
measurement is taken at the same position and just before the 
OFF measurement, and so doesn't involve additional overhead. 
There may be other small overhead losses in moving from row to 
row and starting scans, but these losses can usually be neglected. 

6.2. Integration Time and RMS Noise 
Based on the following definitions: 

77 = Convolution function factor. This is a measure of how much 

the data oversamples the convolution function. 
ti-ow = Time to scan one row in the image. 
Nrow = Number of rows in image. 



8 



J. G. Mangum, D. T. Emerson, and E. W. Greisen: The On The Fly Imaging Technique 



Or„H> = Row separation. 
R = Scanning rate. 

Nceii = Number of independent sampling intervals in image. 
&Nyqiiist = Nyquist sampling interval (= A/(2D)). 
Tsys = System temperature. 
Av = Spectral resolution. 

we can write the relation for the integration time per sampling 
interval as: 



tcell 



row row 

Ncell 
2 

Nyquist 



ROn 



and the RMS noise per sampling interval is: 



(23) 



O'cell - 



ROrow 
TjAv 



(24) 



The factor 77 is dependent upon the convolution function used 
to grid the OTF data. The 77 factors for the convolution functions 
available in the AIPS single dish analysis programs are listed in 
Tabled! 



6.3. Differential Doppier Correction 

Since the correction for the radial motions of the local standard 
of rest relative to a given source vary as a function of position 
on the sky, time, and telescope location, the conversion from 
spectral frequency to velocity varies over a given image. This 
is a subtle effect, traditionally ignored in spectral-line imaging. 
However, for a coordinate 1.4° from the reference (e.g. a 2° 
square field with a center reference), the error can be as large as 
1.16 km s"' (LSR velocities) or 0.79 km s"' (heliocentric) when 
observed with the NRAO 12 Meter Telescope. The LSR velocity 
error changes significantly with time of year. Narrow-band ob- 
servations of wide fields, such as observations of cold molecular 
regions, may be seriously affected by this effect. Fortunately, so 
long as the spectra are fully sampled in frequency, this effect can 
be fully corrected by doing the LSR calculations separately for 
every point of each map, for all dates, and then interpolating to 
a common frame and date in frequency space. 



6.4. Smoothing and Signal-To-Noise for Critically-Sampled 
Data 

When you grid a critically-sampled (Nyquist-sampled) OTF map 
with a sensible convolution function (such as the gaussian- 
tapered Bessel function described in the resulting spa- 
tial resolution is only minimally degraded. However, the noise 
spectrum of your gridded map is flat (or, equivalently, the 
noise in each gridded point is independent), while the astron- 
omy spectrum follows that of the spatial frequency response 
of the antenna. If you now smooth the data significantly - say 
to a resultant beamwidth several times greater than the orig- 
inal spatial resolution - both the astronomical signal and the 
noise get smoothed. The effective increase in integration time is 

given roughly by new convolution function area obviously the 
b b J J origmal samplmg cell area -' 

signal-to-noise ratio increases by the square root of this quan- 
tity. Note that it is the ratio of the new convolution function to 



the original sampling that counts. This is because the high fre- 
quency terms of the original noise, which probably dominate the 
total noise, are reduced in amplitude much more than the astro- 
nomical signal, whose high-frequency components were already 
weighted down. 

After you have smoothed the data once, the noise and astro- 
nomical signal now have a very similar spectrum. If you smooth 
the data a second time, again to a beamwidth several times larger 
than the current (smoothed) beamwidth, the effective integration 

time now increases roughly by nsw beam aiea 

b J J current beam area 
ratio of the beams, independent of the sampling of either mea- 
surement. This in a sense is because the noise in adjacent sam- 
ples of the smoothed data is no longer independent, so you don't 
gain as much in the second smoothing as was gained in the first. 

The conclusion is that to be rigorous you need to know the 
noise spectrum as well as the data beamwidth in order to cal- 
culate what the effective integration time becomes when you 
smooth critically-sampled data. It depends on whether you can 
consider noise in adjacent samples independent. After the initial 
gridding of a critically-sampled data set, the noise is independent 
in adjacent samples. After you have smoothed the data once, the 
noise is no longer independent. 



7. OTF Observing at the NRAO 12 Meter Telescope 

Before starting a spectral line or continuum OTF map at the 12 
Meter, the telescope control system was configured to raster map 
the target field. For spectral line OTF, the map was taken in a 
total power observing mode in the sense that you acquired a cal- 
ibration spectrum (a vane calibration) and a total power off mea- 
surement, followed by one or more total power scanning rows, 
typically made up of hundreds of individual spectra (the on mea- 
surements). For continuum OTF, the map was acquired using 
the continuum beam-switched observing mode. In this mode, the 
subreflector is switched between two positions (the and "-" 
beam) in azimuth while the telescope scans. Each of the indi- 
vidual spectra or continuum "+" and "-" beam total power mea- 
surements is tagged with the actual antenna encoder positions. 
As a result, antenna tracking errors caused by wind gusts, for ex- 
ample, are actually recorded and taken into account in the data 
analysis stage. 

For spectral line OTF, the same calibration and off measure- 
ments were used to calibrate all of the on measurements in the 
scanning rows until another off measurement was taken. Each 
total power off measurement was given its own scan number. All 
of the spectra or continuum "+" and "-" beam total power mea- 
surements in each scanning row were concatenated along with 
arrays of time and position information and stored on disk under 
a single scan number with a single header The header informa- 
tion for each scanning row contained the scan number of the pre- 
vious off (for spectral line OTF) and calibration measurements. 

Given the scanning parameters, the positioning system was 
configured for tracking rates, row offsets, and the duration of a 
scanning row. The spectral line data taking backend instructed 
the tracking system to move to position and begin scanning. 
Using some handshaking bits on the telescope's status and mon- 
itor (SAM) bus, the two systems were synchronized at the start 
of each row. In addition, both the tracking and backend com- 
puters had IRIG clock cards that were driven by the observatory 
GPS standard clock. The data backend read out the spectrometer 
system every 100 ms for spectral line OTF. The continuum data 
system was read out every 250 ms. For both spectral line and 
continuum OTF, the data backend tagged each data parcel with 
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Function 



Farms 



Function Form 
Linear Circular 


R 


1.00 


0.78 


0.5 


3.14 


3.14 


3.0 


1.36 


1.16 


3.0 


2.33 


1.43 


3.0 


3.59 


3.03 


3.0 



Pill Box 
exp[-(My 



Ml) 



X exp 



X exp 



b= 1.00 
a= 1.14 

a= 1.55, b=2.52 
b = 2.52, c = 1.55 




CO 1-0 



a (1950) 

Fig. 13. CO 1 — > integrated intensity image of a molecular 
cloud near the Pleiades cluster The Pleiades member Merope is 
located in the upper left portion of this image, just off of the main 
emission component. The CO emission appears to be tracing an 
interaction between the young stars in the Pleiades cluster and a 
nearby molecular cloud. 



the UT time stamp. For the digital autocorrelation spectrometer 
at the 12 Meter Telescope (called the Millimeter Autocorrelator 
(MAC)), before the data could be read out Fourier transforms of 
the 100 ms data samples needed to be calculated. To make data 
analysis as fast as possible, the FFT's were performed in real- 
time by the MAC control computers. Every 10 ms, the tracking 
computer recorded its Az/El position with respect to the field 
center Finally, at the end of the scanning row, the position infor- 
mation from the tracking system was merged with the data. An 
interpolation of the position information was then made to align 
slight differences between the time stamps of the two data sets. 



8. AlPS Processing 



AIPS ( Greisen| ( 2003 i) was the default analysis software used 
to process OTF data acquired with the NRAO 12m Telescope 
(see Mangum ( |1999 1). The main processing task was SDGRD 
(developed by Eric Greisen), which read the raw 12m ("sdd") 
data and gridded it using a variety of user-selectable convolu- 
tion functions. Figure 13 shows a sample of an OTF image ob- 



tained with the NRAO 12 Meter Telescope and processed with 
the AIPS OTF image analysis tasks. 



9. Conclusion 

The On-The-Fly (OTF) imaging technique is an efficient single- 
dish radio telescope imaging mode that has been implemented 
at several radio telescope observatories. The relative gain in ob- 



serving efficiency is most extreme for the "classic" radio tele- 
scopes still in operation. This was certainly true for the NRAO 
12 Meter Telescope. The NRAO 12 Meter Telescope began life 
as the 36 Foot Telescope, the telescope responsible for the birth 
of millimeter-wavelength molecular astronomy. Its history is one 
of success, innovation, and an unparalleled desire by the staff 
operating this facility over its 32-year lifetime to provide a high 
level of service to the astronomical community. 

The 36 Foot Telescope made its first millimeter-wavelength 
astronomical measurements in October 1968. This was the start 
of a period of explosive growth in this new area of astronomical 
research, during which most of the dozens of molecular species 
known to exist in the interstellar medium were first detected with 
the 36 Foot. 

In the early 1980s the telescope's reflecting surface and sur- 
face support structure were replaced and in 1984 it was re- 
christened as the 12 Meter Telescope. Its scientific program sub- 
sequently evolved from one dominated by astrochemistry to a 
broader mix of studies of molecular clouds and Galactic star for- 
mation, evolved stars, astrochemistry, and external galaxies. 

The NRAO 12 Meter Telescope was the only millimeter- 
wavelength telescope in the U.S. operated full-time as a national 
facility. More than 150 visitors used it each year It offered users 
flexibility and the opportunity to respond quickly to new sci- 
entific developments. Its low-noise receiving systems covered 
a wide range of frequencies - all atmospheric windows from 
68 GHz to 300 GHz - and much attention was given to mak- 
ing the instrument work reliably throughout this range. Flexible 
spectral line and continuum back-ends allowed the observer to 
match the instrument to the scientific goals. The development of 
multi-beam receivers and the on-the-fly observing technique in- 
augurated a new era of high-speed source mapping on angular 
scales complementary to those of the millimeter-wave interfer- 
ometers. The telescope control system offered great flexibility 
and provided a proven remote observing capability. 

On February 22, 2000, NRAO announced that it would close 
the 12 Meter Telescope at the end of the current observing sea- 
son. On 26 July 2000, the NRAO 12 Meter Telescope made its 
final astronomical measurements as a U.S. national astronomi- 
cal facility. The telescope however continues to be operated as 
part of the University of Arizona's Steward Observatory. The fi- 
nal NRAO data set was an on-the-fly image of the CO 1 — > 
emission from a star formation region in the Cygnus-X region. 
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